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Abstract 

We propose a framework for quantitative evaluation of dynamical tendency for polar- 
ization in arbitrary random variable that can be decomposed into a pair of orthogonal 
subspaces. The method uses measures based on comparisons of given dynamics to its 
counterpart with statistically independent components. The formalism of previously 
considered X -distributions is used to express the aforementioned comparisons, in ef- 
fect putting the former approach on solid footing. Our analysis leads to definition of 
a suitable correlation coefficient with clear statistical meaning. We apply the method 
to the dynamics induced by pure-glue lattice QCD in local left-right components of 
overlap Dirac eigenmodes. It is found that, in finite physical volume, there exists a 
non-zero physical scale in the spectrum of eigenvalues such that eigenmodes at smaller 
(fixed) eigenvalues exhibit convex X-distribution (positive correlation), while at larger 
eigenvalues the distribution is concave (negative correlation). This chiral polarization 
scale thus separates a regime where dynamics enhances chirality relative to statistical 
independence from a regime where it suppresses it, and gives an objective definition to 
the notion of "low" and "high" Dirac eigenmode. We propose to investigate whether 
the polarization scale remains non-zero in the infinite volume limit, in which case it 
would represent a new kind of low energy scale in QCD. 



1 Introduction 

In physics it is sometimes useful to decompose a multi-component variable/quantity with 
respect to a pair of orthogonal subspaces spanning its range. Examples relevant for the 
area of application that we target here, namely QCD and its vacuum structure, include the 
decomposition of a bispinor into left-right chirality components, and the decomposition of 
the gauge field strength tensor into self-dual and anti-self-dual parts. Physical reasons for 
considering such decompositions usually have to do with the fact that the dynamics of the 
variable involves (or is expected to involve) a specific type of relation among such suitably 
defined projections. The expected correlations can then be tested in a relevant physical 
experiment or evaluated in theoretical calculations. 
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The property of polarization, namely the tendency for asymmetric participation of the 
two subspaces in the preferred values of the variable, is frequently of interest in the above 
context. While the qualitative meaning of the concept is obvious, the ways to characterize 
and quantify the degree of polarization in the dynamics are highly non-unique a priori. One 
typically operates under the implicit assumption that the physical context uniquely "selects 
the tools", i.e. that the choice of the polarization characteristics is driven by the physical 
nature of the quantity at hand. There are situations however, the study of QCD vacuum 
structure being one of them, where we may be interested in polarization characteristics that 
pertain to dynamics alone and are not conditioned on which specific physical situation the 
dynamics is modelling. Our goal in this paper is to develop a framework to describe and 
quantify such dynamical tendency for polarization. In addition to being free of kinematic 
effects, we also require that the framework provides a possibility for high level of detail. In 
other words, we are not only interested in a suitable average polarization parameter, but 
rather want the aspects of polarization to be separated from other content of the dynamics. 

To construct our formalism, we will build upon the X-distribution proposal of Ref. pQ 
which grew from similar motivations albeit with a specific aim to characterize the local 
chirality of QCD Dirac eigenmodes. The X-distribution approach is based on the following 
logic. Consider the left-right decomposition ip = ipL + ipR of the local value in the eigenmode, 
and some function X = X{^Li^r) that reflects the degree of chiral polarization in ip. Since 
X is defined at every space—time point, it is a local characteristic providing the information 
on changes in the left-right polarization within the mode. The equilibrium ensemble of 
QCD gauge configurations induces the distribution of associated X-values for chosen sets 
of eigenmodesQ In other words, denoting by X an independent variable whose domain is 
in the range of X, there exists a function -P(X), namely the probability distribution of the 
polarization characteristic specified by X. The X-distribution -P(X) provides rather detailed 
information about the chiral polarization properties of the selected group of modes. 

While the concept of X-distributions is suitable for our purposes, the issue of "dynam- 
ical versus kinematical" has not been explicitly addressed in that framework. The channel 
through which kinematics enters here is the function X since there is a large freedom in 
choosing measures that qualitatively reflect chiral polarization. Rather than analyzing this 
freedom, the related work was carried out with the implicit rationale that the proper fixed 
choice of this polarization function can be made. For example, the original definition in 
Ref. [1] involved X linear in the polar angle of point I^rD from the "left-right plane" 

(chiral orientation parameter). This was used in most of the early follow-up works [2 J with 
some departures [21 H] • However, especially in the case of Dirac eigenmodes, it is not entirely 
clear which physical consideration should fix the selection. Thus, focusing exclusively on the 
dynamics, in the sense discussed above, provides a meaningful direction to proceed in. 

With the above discussion in mind, we wish to propose a construct that (1) provides 
a detailed differential information on polarization modeled after the X-distribution but (2) 
doesn't rely on a fixed polarization function X, and (3) reflects only the dynamics. Note that 
while the involvement of kinematics proceeds via the polarization function X, requirement 
(2) only represents a necessary condition for this characteristic to be considered purely 

1 In the early works, additional constraints to define sub-ensembles of X-values, such as restrictions to 
space-time regions with strong fields, were frequently used. 
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dynamical. Indeed, we view the concept of "dynamical" in this context as entailing a specific 
meaning of "correlational" or, in more general terms, based on some form of comparison to 
statistical independence. This is why the requirements (2) and (3) are listed separately and 
need to be fulfilled simultaneously. 

To satisfy these requirements, we first reinterpret the X-distribution as a relative mea- 
sure. In other words, we show that one can view -P(X) not only as an object assigned to given 
dynamics and the polarization function, but also as an object comparing this dynamics to 
another dynamics associated with the polarization function in question. We develop a formal- 
ism that allows us to switch freely between the two interpretations which, in technical terms, 
allows us to construct an X-distribution for arbitrary pair of dynamics. Equipped with this 
tool, we then proceed to define a truly dynamical characteristic as a relative X-distribution 
assigned to dynamics of interest and the related one where left and right components are 
chosen independently of one another from corresponding marginal distributions. Note that 
the latter is the needed point of comparison representing statistical independence. The cru- 
cial point, implicitly embedded in this procedure, is that it is independent of the polarization 
function used to carry it out. In other words, while the implementation must proceed via 
some polarization function X, the final result doesn't depend on this choice. We thus refer 
to the above dynamical characteristic as an absolute X-distribution in order to distinguish 
it from the original one that can only be viewed in a relative sensed 

In terms of requirement (3), one can draw a parallel between our construction and the 
definition of the Pearson correlation coefficient for random variables qi, q2 governed by a joint 
probability distribution V(qi,q2), namely r = (((71(72) — (<7i) (92)) / '(ci ^) , with a±, o~i being 
the corresponding standard deviations. Computing just ((71(72) will not convey any dynamical 
information (it won't tell us anything about tendencies of the relationship between qi and q^)? 
nor does this value have any specific statistical meaning. However, after the subtraction of 
"uncorrelated relationship" (<7i)((72), the resulting covariance acquires a qualitative dynamical 
content. Finally, the normalization via standard deviations, ensuring that — 1 < r < 1, 
gives the coefficient a well-defined statistical meaning, and thus a quantitative dynamical 
significance. In effect, the dynamical construct that we propose implements the "subtraction 
of statistical independence" at the level of distributions which, among other things, has an 
added bonus that proper normalizations are automatically in place. 

With absolute X-distribution at hand, one can naturally define an averaged quantity 
based on it, and this role is played by the correlation coefficient of polarization (CCP). CCP 
has a clear statistical meaning in that it is linearly related to the probability that the sample 
chosen from the distribution governing the dynamics is more polarized than the sample cho- 
sen from its counterpart with statistically independent components. One utility of CCP is 
that it naturally ranks (orders) possible dynamics with respect to the polarization. In partic- 
ular, larger CCP means larger dynamical tendency for polarization. Moreover, the positive 
CCP signals that the dynamics enhances the polarization relative to statistical independence, 
while the negative value means that polarization is being dynamically suppressed. Due to 
its absolute nature, we propose CCP as a basic "figure of merit" quantifying polarization 
properties inherent to dynamics. 

2 One can also refer to this construct as the dynamical X-distribution, explicitly emphasizing the fact that 
its purpose is to remove any kinematical content from the description of polarization. 
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As a first application of the above techniques, we analyze local chiral properties of low- 
lying QCD Dirac eigenmodes. It was already the motivation of Ref. [T] that local properties of 
these modes are expected to be connected to low-energy features of QCD, especially to those 
related to spontaneous chiral symmetry breaking. However, our goals in this investigation are 
quite different from those of Ref. [1] . Rather than checking the consistency of assumptions for 
certain models, as we did in Ref. [5] which contains the preliminary discussion of techniques 
developed here, our aim is to build a model-independent information on relevant properties 
of these modes in line with the bottom-up approach to QCD vacuum structure as formulated 
in Ref. [6]. Indeed, it is the purpose of this series [7] to contribute to developing the set of 
techniques that can facilitate these goals. 

The summary of our main initial findings on the chiral polarization properties of low- 
lying overlap Dirac eigenmodes are as follows. The absolute A-distributions of the lowest 
modes are convex on their domain, while the A-distributions of higher modes are concave. 
This implies that the lowest modes exhibit a dynamical tendency for chirality (positive 
CCP) while the higher modes dynamically suppress it (negative CCP). At fixed volume and 
lattice spacing, there exists a transition point At(V, a) in the spectrum where the absolute 
A-distribution is strictly flat, and the polarization properties are indistinguishable from 
statistical independence. Based on four lattice spacings, our data supports the proposition 
that this scale remains finite in the continuum limit, i.e. lim^o At(V, a) = At{V) > for 
V < oo. The appearance of this sharply-defined new feature in the QCD Dirac spectrum 
could be relevant for understanding the mechanism of spontaneous chiral symmetry breaking. 

The structure of the paper is as follows. We start by discussing an explicit example 
of the absolute A-distribution calculation in Sec. [2j Our goal here is to provide a concise 
introduction to the method with sufficient practical and conceptual information for it to be 
directly applied to the problem of interest here, namely the study of local chiral properties in 
overlap Dirac eigenmodes. This application and corresponding results are then described in 
Sec. |3j In Appendix |A] we describe in detail the theoretical structure, the meaning and the 
scope of absolute polarization methods proposed here. This is a rather extensive topic of its 
own, and the section can be read independently of the QCD content. Indeed, this discussion 
is carried out in a completely generic setting, emphasizing the universal applicability of the 
formalism to any area where the concept of polarization plays relevant role. 

2 Absolute X— Distribution — a Primer 

Given that our first application of absolute polarization methods targets chiral properties 
of Dirac eigenmodes, we will use the specifics of that setting to introduce the approach. 
For definiteness, consider the set of lattice Dirac eigenmodes comprised of lowest two (non- 
zero) conjugate pairs on gauge backgrounds from ensemble E\ of Table [T][f| Collecting all 
local values if) from these modes into single pool, let us consider the associated collection of 
two-component objects ((71,92); namely the magnitudes of left and right components of if) 

(<Zl,?2) = {\^l\Mr\) lfj L = \{l- l5 )lf) lP R =\{l+ lb )lP (1) 

3 Details concerning the parameters of lattice gauge ensembles and lattice Dirac operator are given in 
Sec. [3] and are not important for this discussion. 
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Figure 1: Left: Dynamics governing the polarization properties of two lowest modes for 
ensemble E\. Right: The associated "uncorrelated" dynamics; solid gray lines separate the 
quadrant into 10 sectors, each containing 10% of the population. 

The distribution of these "samples" is shown on the left side of Fig. [TJ and statistically 
approximates the underlying base probability distribution P&((7i, (72) describing the limit of 
infinite gauge statistics. 

One can assess a degree of polarization in a given sample by its polar angle ip G [0, 7r/2] 
which is convenient to symmetrize into reference polarization coordinate, namely [TJ 

4 

x = -ip - 1 x G [-1,1] (2) 

7T 

Note that strictly left-polarized and strictly right-polarized samples are characterized by 
values x = — 1 and x = +1 respectively. The probability distribution P r (x) of polarization 
coordinate over the population of samples is sometimes invoked as a detailed measure char- 
acterizing chiral polarization properties of Dirac modes in question, and is usually referred 
to as the (reference) X-distribution [TJ. Our goal is to define a construct of similar nature 
which however focuses on the relation of the dynamics in question to the "uncorrelated" one 
where left and right components of if) = (if)L, ipn) are paired with one another at random from 
available pool of possibilities. Indeed, such comparative polarization measure could then be 
viewed as dynamical. Note that the needed point of comparison (the population of statis- 
tically independent components) associated with our exemplary dynamics is shown on the 
right side of Fig. [TJ It approximates the underlying distribution V^(qi, (72) = p(qi)p(q2) with 
p(q) being the distribution of a single component in original Vb(qi,q2)- This uncorrelated 
dynamics carries its own reference X-distribution P r M (x). 

To obtain the polarization measure of the above kind, we perform a differential compar- 
ison of polarization in Vb(qi, (72) relative to P^((7i, (72) as follows. A ray passing through the 
origin can be specified by its reference polarization coordinate x, as shown in Fig. [TJ Deter- 
mine the fraction of population contained between the gi-axis and this ray, separately for 
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Figure 2: Steps in the construction of absolute X-distribution for data shown in Fig. [TJ See 
the discussion in the text. 

Vh and V%. These fractions represent the cumulative probability functions S r (x) and S^(x) 
associated with P r (x) and P^(x) respectively. Eliminating x, every line is represented by a 
point in the cumulative distribution plane (S", S r ) with the set of all such points forming 
a graph starting at (0,0) and ending at (1,1). If the the two situations ("dynamics") in 
question had identical polarizations, the graph would be a straight line. On the left side of 
Fig. [2] we plot the result of this construction for data displayed in Fig. [TJ One can see that 
the correlated and uncorrelated dynamics have very similar polarizations since the graph is 
close to being linear. To obtain a differential comparison, and to better see the differences, 
we compute the slope of the cumulative polarization graph and show it on the right side 
of Fig. [2j Since this dependence represents a probability distribution in variable S", we 
can see that the correlated case exhibits a small excess of probability with respect to the 
uncorrelated one (the horizontal line) near the extremal values. Noting that the extremal 
values SI? = 0, 1 actually correspond to strictly polarized CclSCS X — 1, 1 respectively, the 
observed qualitative behavior of the slope graph in fact conveys the excess of polarization in 
dynamics described by Vb{Qi,Q2) relative to dynamics described by V^{(li,<h)- To standard- 
ize this polarization characteristic, we finally convert the distribution in into distribution 
in rescaled polarization-like variable — > 2S™ — 1 = X whose domain is [—1, 1]. The 
resulting construct Pa{X) will be referred to as the absolute X-distribution. 

We now wish to make several remarks regarding the above method. 

(i) Performing a practical calculation of absolute X-distribution with uniform resolution in 
X requires that we evaluate cumulative probabilities at values Xi that split the uncorrelated 
population into Nj, bins of equal size. This is indicated on the right hand side of Fig. [TJ for 
Nf, = 10 by the rays enclosing 10% of population each, and similarly by vertical lines in 
Fig. [2J The corresponding reference distribution P r (X = x) and absolute X-distribution 
Pa(X) computed with this resolution are shown in Fig. [3J 

(ii) As one can readily inspect, the above geometric construction would lead to the same re- 
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Figure 3: Reference X-distribution (left) and absolute X-distribution (right) for data shown 
in Fig. [I] with resolution of 10 bins. 



suit had we decided to quantify the degree of chiral polarization in ip by some arbitrary polar- 
ization function X(x) rather than by the default choice X(x) = x. Indeed, such different pos- 
sibilities correspond to different parametric representations of the same curve in the S r ) 
plane, and this information drops out of the calculation. The associated reparametrization 
invariance is significant since it expresses the fact that the proposed description of polariza- 
tion is independent of the "language" (or "reference frame" ) chosen to represent the concept 
at the sample level. This is necessary for the description to be considered purely dynamical 
and reflects the absolute nature of this construction. 

(Hi) The definition of absolute X-distribution given above is equivalent to constructing the 
X-distribution of Ref. [TJ, but associated with polarization function 

X A (x) = 2 j* d V F?(y)-l (3) 

This makes the fact that the absolute X-distribution represents a valid polarization mea- 
sure more explicit. In fact, depending on the dynamics Vb(qi,q2) in question, the proce- 
dure adjusts the polarization characteristic so that it directly measures deviations relative 
to uncorrelated dynamics. In this language, the reparametrization invariance of absolute 
X-distribution follows from the fact that is invariant under the choice of reference po- 
larization coordinate as one can directly verify. 

(iv) We wish to emphasize that the strong dependence on the choice of the polarization 
function (dependence on "kinematics") provided a major motivation for extending the old 
X-distribution approach in the direction described here. Indeed, as discussed extensively 
in Appendix |XJ even qualitative conclusions can be entirely misleading if one relies only on 
the X-distribution associated with fixed polarization function. This is illustrated in Fig. [TT1 
where a suitable choice of polarization function can either produce seemingly clear chiral 
double-peaking or what appears as completely non-chiral single peak. 

(v) In addition to detailed description of polarization contained in absolute X-distribution, 
we also define a single figure of merit reflecting the overall dynamical tendency for polariza- 
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tion, namely 

C A = 2 J dX \X\ P A {X) - 1 (4) 

The positive value of this correlation coefficient of polarization Ca £ [— 1, 1] indicates that 
the dynamics enhances polarization relative to randomness, while the negative value implies 
its dynamical suppression. Indeed, as one can check (see Appendix [AJ, Ca is linearly related 
to the probability that the sample chosen from Vb(qi, 92) is more polarized than the sample 
(independently) chosen from V^(qi, 92), leading to the above interpretation. 

The elements described in this section are at the core of the general dynamical polar- 
ization method we are proposing in this work. Formal definitions and detailed discussion of 
the logical structure involved are given in Appendix |A] While the concise overview given in 
this section is sufficient to apply the method to the case of QCD Dirac eigenmodes in what 
follows, we emphasize that the content of Appendix |A] is important not only for deeper un- 
derstanding the scope and the meaning of the approach, but also for its practical application 
in case of singular polarization dynamics. 

3 QCD Dirac Eigenmodes 

We will now apply the absolute polarization methods to the case of low-lying QCD Dirac 
eigenmodes in a more systematic manner. As argued extensively e.g. in Ref. pQ, the proper- 
ties of these modes are expected to encode many important features of QCD vacuum. Their 
direct connection to spontaneous chiral symmetry breaking, beyond the Banks-Casher rela- 
tion [Hj, is particularly anticipated. 

We perform our exploratory calculations in the context of pure glue lattice QCD with 
Iwasaki gauge action [9] . The left side of Table CD summarizes the parameters of our five 
ensembles. Lattice scale has been determined from the string tension following the methods 
and results of Ref. [10J. Ensembles E1-E4 have fixed physical volume, and are intended for 
investigation of the continuum limit. Ensemble E§ has larger physical volume and serves 
mainly to check whether our qualitative conclusions remain stable in that regard. 



Ensemble Size A config Volume Lattice Spacing Aj$$ A£ 



Ei 


8 4 


100 


(1.32 fm) 4 
(1.32 fm) 4 
(1.32 fm) 4 
(1.32 fm) 4 
(1.76 fm) 4 


0.165 fm 


449 


226 


1956 


1980 


E-2 


12 4 


97 


0.110 fm 


407 


169 


1711 


1735 


E3 


16 4 


99 


0.0825 fm 


304 


142 


1513 


1553 


£4 


24 4 


96 


0.055 fm 


344 


136 


1338 


1366 


E 5 


16 4 


99 


0.110 fm 


162 


58 


1087 


1123 



Table 1: The summary of five ensembles used in overlap eigenmode calculations. The right 
side of the table describes some properties of the spectra (in MeV) with A^qw denoting the 
average magnitude of lowest near-zero eigenvalue over the ensemble and A^f H denoting the 
same for highest eigenvalue. A^ 1 ^ is the magnitude of the maximal lowest eigenvalue, and 
A.HIGH the magnitude of the minimal highest eigenvalue. 
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Low-lying eigenmodes of the overlap Dirac operator [TT] were calculated on gauge back- 
grounds from the above ensembles. Specifically, we used the massless overlap operator with 
Wilson kernel, the negative mass parameter p = 26/19 (k = 0.19), and Wilson parameter 
r — 1. The eigenmodes and eigenvalues were computed using our own implementation of 
the implicitly restarted Arnoldi algorithm with deflation [12J. The zero modes and about 50 
pairs of lowest near-zero modes □ with complex-conjugated eigenvalues were obtained for 
each configuration. The right side of Table [1] compiles the information on spectral bound- 
aries of computed regions for each ensemble (see the caption). Note that since zero modes 
are exactly chiral, it is the near-zero modes that are of interest here. They are in fact the 
relevant eigenmodes to consider in the context of spontaneous chiral symmetry breaking. 

3.1 Reference and Absolute X— Distributions 

In order to get the basic idea about the dynamical tendency for chirality in the obtained 
eigenmodes, we first focus on the lowest non-zero and highest computed modes. For this 
purpose, we collect into the "lowest" eigenmode group the two lowest near-zero modes for 
each configuration from the given gauge ensemble. Since the eigenmodes of conjugate pairs 
have identical chiral properties, we actually include one representative from each of the two 
lowest conjugate pairs. Similarly, the "highest" eigenmode group contains eigenmodes from 
the highest two pairs. The collection of all local values of (91,92) (magnitudes of left and 
right components) in the grouped eigenmodes represent the dynamics ^(91,92) that we are 
investigating. 

On the left side of Fig. H]we show the reference A-distributions, or reference polarization 
dynamics, for "lowest" and "highest" eigenmodes in ensembles E\—E±. On the right side 
of the same figure are the associated A-distributions for the dynamics with statistically 
independent left-right components. A detailed description of the method used to obtain this 
"randomized" or "uncorrelated" dynamics is given in Appendix [Bj as are the other technical 
details related to the construction of absolute A-distributions. 

Looking at the left side of Fig. HJ the most obvious feature is the narrowing of the 
distributions going from low to higher modes. This is as expected from many studies of 
standard A-distributions even though we should point out that, contrary to most of such 
studies, here we are not making any restrictions based on magnitudes of ip(x). Rather, all 
local values are part of the dynamics, as they should be. 

An interesting aspect of Fig. H] is the comparison of left and right columns, i.e. the 
comparison of pure kinematics (right side) to kinematics combined with dynamics (left side). 
As one can see, the two are barely distinguishable implying that the qualitative appearance 
of reference A-distributions is almost entirely driven by kinematics. In order to expose the 
dynamical polarization tendencies, we construct the corresponding absolute A-distributions 
that effectively remove the features due to kinematics. The resulting distributions are shown 
in Fig. [5] with lowest modes on the left and highest modes on the right. 

One can easily interpret the absolute A-distributions since the uncorrelated dynamics 
is represented by the uniform distribution shown as a solid horizontal line. Thus, for the 

4 The exact number of computed near-zero mode pairs depends on the number of zero modes (topological 
charge) of a given configuration. If No is the number of zero modes then there were 55 — N$ pairs computed. 
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lowest modes, there is an enhancement of samples with largest polarization relative to statis- 
tical independence, and the corresponding suppression of samples with low polarization. In 
fact, one can read from the graph complete differential information on how much enhance- 
ment/suppression there is for arbitrary segment of uncorrelated population. The correlation 
coefficient of polarization (CCP) is obviously positive for the lowest modes, and the dynamics 
governing them exhibits the dynamical tendency for polarization. As can be inspected from 
the right column, this situation reverses for the highest computed modes in all lattice gauge 
ensembles. Indeed, in that case the dynamics suppresses samples with largest polarization 
and CCP is negative. 

3.2 Convexity of Absolute X— Distributions 

A characteristic feature of absolute X-distributions shown in Fig. [5j as well as of all the other 
ones that we computed from our ensembles, is that they are in some sense the simplest they 
can be. For example, while there doesn't seem to be a reason a priori excluding multiple 
intersections of Pa(X) with uniform distribution on each "chiral half" (interval [—1,0] and 
[0, 1]), we have generically seen just one. Note that these crossings define segments in the 
sample space with dynamically enhanced/suppressed chirality. Moreover, consistently with 
the above, the absolute X-distributions for QCD eigenmodes seem to be either purely convex 
or purely concave on their domain [—1,1]. This implies that they are also strictly monotonic 
on each chiral half. 

To give a well-defined meaning to these kinds of observations, one should specify more 
precisely the dynamics of chiral pairs in question. In other words, one should specify how 
are the associated populations defined in terms of selected Dirac eigenvectors so that the 
corresponding characteristics are both physically meaningful and have proper continuum 
limits. In that regard, things simplify for lattice QCD in the infinite volume since one can 
assign an X-distribution to an individual eigenmode in that case, at least for equilibrium 
gauge backgrounds. Moreover, under standard assumptions, the local properties of such 
modes will not depend on the specific equilibrium configuration in question. Thus, in the 
infinite volume, the eigenvalue A is the only label distinguishing the local behavior of Dirac 
eigenmodes in the underlying gauge theory. It is a common practice that the above is also 
mimicked in the finite volume, except that the gauge ensemble average has to be performed 
in that case. In other words, one defines properties of eigenmodes at scale A by inserting 
5 (A — |A|) into the expression for ensemble average of an observable involving a spectral sum 
over individual eigenmodes 

Here we will follow the above logic and all the conclusions regarding the polarization 
properties of Dirac eigenmodes will implicitly refer to dynamics of modes at fixed A in the 
above senseH With that in mind we suggest the following for further investigation. 

Proposition 1 (Convexity): The absolute X-distributions of low-lying overlap Dirac modes 
in SU(3) pure glue lattice gauge theory are strictly convex, strictly concave or constant. 

5 Note that the practical implementation of this, i.e. the selection of appropriate eigenmodes with only 
finite gauge ensembles at our disposal, can vary depending on how much data is actually available. For 
example, our selection of "lowest" modes in the previous section could be viewed as a procedure to determine 
the behavior of eigenmodes at A = in the continuum limit. 
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Figure 4: Reference X-distributions for "lowest" and "highest" modes in ensembles E\-E^ 
are shown in the left column. The associated distributions with statistically independent 
left-right components ( "uncorrelated" ) are shown in the right column. 
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Figure 5: Absolute X-distributions for lowest (left column) and highest (right column) two 
pairs of modes for ensembles E1-E4. Lattice spacing is decreasing from top to down. 
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Figure 6: Lattice spacing dependence of absolute X-distributions for overlap Dirac eigen- 
modes at A = 1200 MeV (top) and A = 420 MeV (bottom). The corresponding correlation 
coefficients of polarization are shown on the right. 



It is naturally expected that this behavior of absolute X-distributions will not depend on 
the choice of the Dirac operator with exact lattice chiral symmetry. We would also like to 
emphasize that the weaker statement, namely that the absolute X-distributions are strictly 
monotonic or constant on each chiral half, might hold not only for low-lying modes but for 
all Dirac modes. 

3.3 Continuum Limit of Distributions at Fixed A 

Assuming that the dynamics of Dirac eigenmodes at fixed A can be assigned a well-defined 
physical meaning, their dynamical polarization characteristics should have a proper contin- 
uum limit. One should point out though that whether this is true or not constitutes the 
test for the notion of eigenmodes at fixed A rather than the test for applicability of absolute 
X-distributions. 

Since our statistics are rather limited, we define the set of modes specifying the dynamics 
at A > in the following way. For each configuration from given gauge ensemble we take 
two modes ipx^ such that |Aj| < A < |Aj+i|. Here we have ordered the non-zero 

low-lying modes by magnitude as Ai, A*, A 2 , A£, • • •• Since in practice one doesn't encounter 
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any degeneracy of non-zero modes except for conjugate pairs, this definition uniquely picks 
out the two modes from the spectrum that are "hugging" the value A assuming that A is 
sufficiently large (A > A%$# from Tabled]). 

In Fig. [6] we show the convergence of absolute X-distributions for eigenmodes at 1200 
MeV (top) and 420 MeV (bottom). Note that in case of the latter, the result for ensemble 
Ei is not applicable since 420 MeV is less than A^qw ■ This is not a problem because 
Ei appears to be too coarse to be used in continuum extrapolations while our aim is to 
provide the example of the very low A. As one can see from plots on the left, the data quite 
clearly suggests point-wise convergence of absolute X-distributions as the continuum limit 
is approached. To see this in a more coarse-grained manner, we show in the right column 
also the lattice spacing dependence of the moment of the distribution (correlation coefficient) 
in both cases. Their qualitative behavior suggests convergence as well. The lines shown are 
fits to the form c\ + Cia A to the points corresponding to E 2 -E^. While the power four is 
suggested by the data at 1200 MeV, the fit should only be viewed to guide the eye. Based 
on the above, we propose the following for further examination. 

Proposition 2 (Continuum Limit): The absolute X -distributions for overlap Dirac modes at 
scale A have a continuum limit in SU(3) pure glue lattice gauge theory. 

3.4 Chiral Polarization Transition 

The comparison of left and right columns in Fig. El as well as the comparison of top and 
bottom parts in Fig. [6] suggest a high degree of regularity in the dynamical polarization 
properties of the low-lying modes. In particular, the dynamics of lowest-lying modes always 
seems to enhance chiral polarization while the dynamics of higher modes tends to suppress 
it. To check if this is indeed a monotonic feature with respect to scale A of the modes, we 
show in Fig. [7] how the absolute V-distribution for gauge ensemble E4 changes from low 
values upwards. As one can see, the transition from convex to concave indeed proceeds in a 
monotonic manner. 

The above behavior is qualitatively similar for all our ensembles and we are thus led to 
conclude that there exists a sharply defined transition point At in the spectrum of overlap 
Dirac eigenmodes, where dynamical polarization properties of these modes qualitatively 
change. The QCD-inherited dynamics governing the lower part of the spectrum supports 
chiral polarization while in the upper part it produces chiral anti-polarization. Apart from 
the mere existence of this dynamical feature, it is quite remarkable that it is not only the 
CCP that changes sign at At, but rather the whole polarization dynamics (absolute X— 
distribution) appears to be indistinguishable (uniform) from statistical independence. This 
suggests that the existence of this chiral polarization transition represents a robust dynamical 
property of QCD Dirac eigenmodes. 

In order to determine A^ for our ensembles, we scanned the behavior of CCP with respect 
to A and looked for the behavior in the vicinity of the point where it changes sign. As shown 
in Fig. El for gauge ensemble E4, the data suggests linearity in the corresponding region. The 
transition point was thus determined from the linear fit in this neighborhood. More details 
about this procedure can be found in Appendix [Bl 
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Having performed the determination of At for ensembles E\—E±, we show in Fig. [9] the 
dependence of on the lattice spacing. We emphasize that the ensembles in question have 
fixed physical volume. Our main interest here is to examine whether Ay(V, a) has a finite 
continuum limit at fixed V . In other words, whether there is a well-defined finite scale 
associated polarization dynamics of QCD Dirac eigenmodes in finite volume. As one can see 
from the figure, our data clearly suggests that this is indeed the case. The fit to C\ + 02a 4 , 
with the coarsest ensemble E 1 excluded from the fit, was again used to guide the eye. Based 
on our results, we propose the following for further investigation. 

Proposition 3 (Chiral Polarization Scale): Consider pure glue SU(3) lattice gauge theory in 
finite physical volume V. There exists a scale At{V, a) in the spectrum of overlap Dirac 
operator such that its eigenmodes are chirally polarized below At(V, a) and chirally anti- 
polarized above it. Moreover, the absolute X -distribution changes from convex to concave 
and is thus uniform at the transition point. The continuum limit lim a _^o At{V, a) exists and 
is non-zero. 

It should be pointed out that the volume of our ensembles is rather small and it is not clear 
at this point whether the above chiral polarization scale remains finite also in the infinite 
volume limit. We consider this an interesting question that should be studied in the future. 

4 Discussion and Conclusions 

Polarization tendencies, whether observed empirically in natural phenomena or studied the- 
oretically in various model systems, are frequently of interest in physics. In fact, they are 
arguably the most basic attributes of dynamical behavior one can attempt to describe. For 
this reason they are also of interest for the bottom-up program of studying QCD vacuum 
structure [5]. Indeed, the major aspect of this approach is to systematically characterize the 
properties of the vacuum in a model-independent manner. 

Polarization is a highly "coarse-grained" concept and, as such, offers a variety of de- 
scriptions that may be consistent with its intuitive meaning. However, such plurality is not 
necessarily a sign of redundancy. Rather, various descriptions may reflect large freedom in 
possible aspects of polarization that one could be interested in exploring. Once the intent 
of the description is meaningfully specified, suitable quantifiers can emerge in a reasonably 
unique manner. In this work we discussed polarization properties from a different angle. 
Indeed, rather than looking for description tailored for a specific physical context, our goal 
was to construct a differential characteristic that is meaningful regardless of what the dy- 
namical variable in question represents. As such, it could then be viewed as an attribute of 
the dynamics alone, a dynamical polarization characteristic. 

To achieve the above, we considered a large set of descriptors consistent with the qualita- 
tive notion of polarization, namely all possible X-distributions labeled by various polariza- 
tion functions. If one imagines a rule of assignment associating a suitable description with 
any given physical context, then the set of X-distributions would be the range of such map: a 
proper description is always found somewhere within this set. However, given our goals, one 
should not view different polarization functions as representing all possible physical contexts, 
but rather as different "reference frames" in which any given dynamics can be looked at in 
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Figure 9: Dependence of chiral polarization point on the lattice spacing at fixed volume. 



terms of polarization. The X-distribution associated with such a frame is kinematical in the 
same sense as velocity v in given inertial frame is kinematical for mechanical particle: it can 
be made arbitrary by choosing the reference frame appropriately. The necessary condition 
for a quantitative measure to be dynamical is its invariance under the change of reference 
frame, as exemplified by dv/dt and its Galilean invariance in case of Newtonian particle. 
One of our main points in this work is that, for situations involving polarization, the analog 
of acceleration is the absolute X-distribution with its reparametrization invariance. 

Similarly to the case of mechanical particle where one could consider e.g. higher time 
derivatives, it is possible to construct other X-distributions with reparametrization invari- 
ance. Thus, frame independence is not sufficient to fix the dynamical polarization charac- 
teristic uniquely. The needed additional requirement we used was that the characteristic be 
"correlational" . Indeed, as emphasized throughout this article, the absolute X-distribution 
represents the differential comparison of polarization in given dynamics relative to the case 
of statistically independent components. If one performs this comparison on average rather 
than differentially, the correlation coefficient of polarization is obtained. This coefficient 
can thus be viewed as descending from absolute X-distribution via a coarse-graining proce- 
dure. It has a well-defined statistical meaning and provides a suitable measure for ranking 
of possible dynamics in terms of their dynamical polarization properties. 

Somewhat different perspective on our proposal is obtained when one views it as a reduc- 
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Figure 10: Absolute (dynamical) and uncorrelated (kinematical) X-distributions for gauge 
ensemble E$ using the reference polarization coordinate for the latter. Results for lowest 
two non-zero pairs are shown on the left while the highest two pairs are shown on the right. 



tion process wherein the full information on orthogonally decomposed variable Q = (Q\, Q2), 
described by the density function Vf(Q), is being restricted into the information relevant for 
describing its polarization properties only. Integrating out unrelated information naturally 
leads to the notion of polarization dynamics or, equivalently, an X-distribution -P(X), but 
the associated non-uniqueness issues lead to the conclusion that this reduction is in fact 
insufficient. One should rather consider the reduction to a pair 



V f (Q) — > V b ( qi ,q 2 ) — > P(X) , P U (X) (5) 

where P U (X) descends from Vf(Q), namely the distribution of statistically independent 
components Q\ and Q 2 - While both members of the above pair vary with the choice of the 
polarization coordinate, there is an invariant information in the "difference" of the two, and 
we thus propose that the polarization should rather be characterized by 



V f (Q) — ► P A (X),P U (X) (6) 

where the first part is invariant (dynamical) and the second part is variant (kinematical) J§ 
Note that the kinematical part is still a useful ingredient in the description of polarization. 
Indeed, if one stores the information on Vf(Q) in the above form with standard choice of the 
polarization coordinate, such as Pa{X) , P"(x) for reference polarization coordinate, then it 
is possible to recover all other X-distributions one might need in various physical contexts 
that this dynamics could be modelling. As an example of such practice, we show in Fig. [10] 
the polarization characteristics for gauge ensemble E 5 and the group of lowest two (left) and 
highest two (right) non-zero modes computed. The above conceptual aspects as well as the 
theoretical structure underlying the method are discussed at length in Appendix [A] 



6 Note that, as opposed to pair in (O, the arguments of these two ^-distributions have different meanings 
in this case. 
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Taking a QCD perspective, it is interesting to apply the above techniques to Dirac eigen- 
modes and investigate their local chiral structure. This is where the related issues first 
arose jl]. To this end, we obtained and presented the basic picture of dynamical polarization 
properties in low-lying overlap Dirac eigenmodes for pure glue QCD. The most intriguing 
aspect of our findings relates to emergence of new sharply defined feature in the Dirac spec- 
trum, namely that of a chiral polarization transition. Indeed, our data indicate that, in finite 
physical volume, there exists a scale in the Dirac spectrum such that QCD dynamics leads 
to enhancement of local chirality in modes at lower scales, while the same dynamics works to 
suppress it in modes above the transition. In substantiating this, we focused on continuum 
limit at fixed physical volume but this basic feature appears to be robust when the volume 
is increased, as indicated by the results for ensemble E 5 shown in Fig. [TU1 

Even though it has long been suspected that chirality properties of low-lying Dirac modes 
are relevant in the mechanism of spontaneous chiral symmetry breaking, it was not known 
before that QCD dynamics leads to such a clear-cut qualitative distinction in the behavior of 
low and high parts of the Dirac spectrum. The above result not only confirms the conceptual 
value of the approach to polarization that we developed here, but also offers what could be 
a fruitful avenue for "bottom-up" investigation of spontaneous chiral symmetry breaking in 
QCD. An immediate point of interest in this regard is whether the chiral polarization scale 
At remains finite in the infinite volume. Indeed, the existence of such chirality-related scale, 
generated by QCD dynamics, would open an interesting new class of possibilities for the 
mechanism of spontaneous chiral symmetry breaking. 
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A General Description of the Method 

In this section we describe the absolute X-distribution method in detail. Since the formalism 
is applicable in a fairly generic setting, we carry this discussion out without reference to a 
specific physical situation in question. This helps to bring the essential ingredients more 
to the forefront. Consider a quantity Q taking values in linear space with scalar product. 
Assume further that there is some physically motivated decomposition of Q into a fixed pair 
of orthogonal subspaces in place, namely 

Q = Qi + Q 2 Qi • Q2 = (7) 

where " ■ " denotes the scalar product. What we implicitly have in mind here is that the or- 
thogonal subspaces in question are "equivalent" e.g. of the same dimension, but the resulting 
methods can be generalized to include non-symmetric cases. 
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What we mean by specifying the "dynamics" of Q is selecting an arbitrary probability 
space with allowed values of Q as its sample space. For practical considerations (and ma- 
nipulations) we will assume however the existence of an underlying probability distribution 
(function or generalized function) Vf(Q) from which the probability of measurable events can 
be obtained via corresponding integrations. It is not important for this discussion whether 
and how this distribution descended from dynamics of a more complex or fundamental sys- 
tem. For example, in applications related to studying space-time structure of composite 
fields in QCD (QCD vacuum structure) this probability distribution would have ultimately 
descended from probability density of gauge configurations given by QCD action. 

Since our goal is to describe polarization properties with respect to (EJ), it is convenient to 
consider dynamics in terms of the polarization components, i.e. Vf(Qi,Q 2 ). In addition to 
dynamics, the method we will be describing is built on the notion of a polarization function 
'X(Q 1 ,Q 2 ). The role of X is to quantify the degree of polarization for any sample (Qi,Q 2 ). 
One can think of it in terms of a "polarization meter" accepting a sample as an input and 
providing its polarization value as an output. Different devices of this type are characterized 
by different polarization functions. In what follows, we will carry out our discussion in the 
situation where the dynamics are symmetric, namely Vf(Qi,Q 2 ) = Vf(Q 2 ,Qi), and the 
polarization meters are also symmetric, namely |X(Qi, Q2) | = |X(Q2, Qi)\- This is most 
relevant for our targeted applications. However, we will keep the framework and notation 
wide enough to incorporate future generalizations in this direction. 

Finally, we will be quantifying the tendency for polarization by relative magnitudes of Q\ 
and Q 2 . For example, by saying that Q is strongly polarized in the first direction we mean 
that \Qi\ 3> |6?2|- Consequently, the probability distribution of component magnitudes, 
namely 

V b ( qi ,q 2 ) = j dR 1 dR 2 V f (R 1 ,R 2 )5(q 1 - ^1) 5(q 2 - \R 2 \) q t = \Qi\ (8) 

will be the object characterizing "dynamics" for our purposes, with the subscript "6" in Vb 
indicating that this is the representation of dynamics in base sample space variables (qi, q 2 ). 

A.l Polarization Functions 

To start with the X-distribution method, one needs to select a suitable function X(Qi,Q 2 ) 
serving as a measure of polarization on the "sample space" . Such polarization function will 
thus reflect the relative contribution of the two orthogonal subspaces in Q. Due to the linear 
structure involved, the values of X assigned to Q and its arbitrary rescaling ctQ [a 7^ 0) 
should be identical, i.e. we consider Q and aQ to be equally polarized samples. Restricting 
ourselves to characteristics X(qi,q 2 ) that only depend on magnitudes (91,52) this translates 
into the requirement that X(q\, q 2 ) only depends on t = qi/qi- 

Note that, in terms of t, "being smaller" indicates being more polarized with respect 
to the first direction while "being larger" signals larger polarization with respect to the 
second direction. Thus, even a simple choice X(qi,q 2 ) = q 2 /qi has some basic features we 
are seeking. However, we wish to work with characteristics that are symmetric with respect 
to the two subspaces, and have a standardized compact range. To achieve this, we impose 
the following set of requirements on acceptable functions X. 
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Definition 1 (Polarization Functions): Consider real-valued functions X(q\,q2) on M.q x 
where M.^ = [0, oo). The subset of such functions satisfying the following conditions: 

(a) X(q,tq) is non-decreasing function of t G M.q for fixed q > 

(b) X(0,g) = 1 for q > 

(c) X(q,tq) is independent of q > for fixed t G Mq 

(d) X(q 1 ,q 2 ) = -X(q 2 ,qi) 

will be referred to as a set of polarization functions. 

Note that, with any polarization function, samples (51,92) strictly polarized in the first 
direction (qi > 0,q 2 = 0), strictly unpolarized (gi = q 2 ), and strictly polarized in the second 
direction (g x = 0,q 2 > 0) are characterized by functional values —1, and +1 respectively. 
X serves as a comparator in a sense that we say that sample (qi, q 2 ) is more polarized than 
sample (r 1; r 2 ) if \X(qi, q 2 )\ > |X(ri, r 2 )\. Similarly, given our conventions, we consider sample 
(gi, q 2 ) more polarized in the first direction than sample (r 1; r 2 ) if X(qi, q 2 ) < X(r 1; r 2 ) while 
we consider it more polarized in the second direction if X(gi,g2) > X(ri,r 2 ). There are 
obviously infinitely many such comparators (polarization functions). 

Every polarization function vanishes at origin (0, 0) and is otherwise specified by the 
choice of non-decreasing X(q 2 /q\) = X(t) such that X(t) = — X(l/t) and X(0) = —1. Variable 
t is in one to one correspondence with the polar angle tp = tan _1 (t) in Cartesian plane, and 
it can thus be useful to parametrize x R^" using polar coordinates (p, y?). In this case 
condition (c) translates into X(p, tp) = X(ip). Further simplifications arise by "symmetrizing" 
the angular variable tp G [0,7r/2] via rescaled polar angle x of Eq. (T5]). In this case the 
remaining conditions on non-decreasing X(x) translate into X(x) = —X(—x) and X(l) = 1. 
In other words, the set of polarization functions can be faithfully parametrized by odd, 
non-decreasing functions on [—1, 1] with unit maximal value. 

When viewed as functions of x, then simple examples of polarization functions are odd 
integer powers and their generalizations to positive real powers, namely 

X c {x;a) = sgn(x) \x\ a a>0 (9) 

Note that X c (x; 1) = x corresponds to the original chiral orientation parameter of Ref. [TJ. 
When viewed as functions of t, then another simple family of polarization functions is spec- 
ified by 

X R {t;a) = a > (10) 

t ~\~ 1 

The special case of X R (t; 2) was considered in [3]. Another possibility is to work with 

X G {t]a) = -t&n-\t a ) - 1 a>0 (11) 

7T 

where a = 1 is the chiral orientation parameter of Ref. [T] and a = 2 was used in [I]. 
A.l.l Set X 

Let us denote the set of all odd non-decreasing real-valued functions on [—1, 1] with unit 
maximal value as Xq. Each polarization function is labeled by an element of Xq and, in 
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fact, we will write X G Xo regardless of whether we think of X in terms of a single variable 
), or in terms of two variables X = X(q 1 ,q 2 ). 
Before proceeding, it is useful to make a subtle distinction that will be needed later. Since 
functions in Xo are bounded and non-decreasing, they are continuous except possibly on a 
countable set of points where they have finite jumps. We will group the functions that only 
differ by values at points of discontinuity and treat them as a single measure for polarization. 
To describe the resulting equivalence classes as a new function set, we represent each class 
by a function that is middle- valued for \x\ < 10 We then have the following representation 
of polarization functions. 

Definition 2 (Set X ): Set X is a subset of real-valued functions f(x) on [—1, 1] satisfying 
the following conditions: 

(1) f(x) is non-decreasing 

(2) /(I) = 1 

(3) f(x) is middle-valued at \x\ < 1 

(4) f(-x) = -f{x) 

It should be emphasized that even though X C X , it implies a very mild restriction on the 
set of polarization functions. Indeed, for every element of X there is an element in X which 
differs from it at most on a countable set in domain [—1, 1]. The rationale for this restriction 
is that polarization functions represented by elements of X are in one to one correspondence 
with possible polarization dynamics as will be discussed shortly. 

A. 2 Polarization Coordinates and X— Distributions 

When studying the polarization properties of dynamics Vb(qi, q-z) using polarization function 
X(qi,q 2 ), it is convenient to express probability dependencies directly in terms of X. In 
particular, we define the following associated objects 

V(X,p) = J™ dq t jT dq 2 V b {qi, q 2 ) d(x - X{q± , g 2 )) s(p - yjql + gf ) (12) 

P{X) = J d Ql ^ dq 2 V b (q 1 ,q 2 )5( y X-X(q ll q 2 ^ (13) 

poo POO 

r = / dq 1 dq 2 V b {qi,q 2 ) | X(q 1 ,q 2 ) | (14) 
Jo Jo 

where X denotes a generic independent variable parametrizing the range of polarization 
functions, i.e. X G [—1,1], and p G is the magnitude of sample Q. In this way, every 
polarization function X(qi,q 2 ) G X defines a corresponding "sample space coordinate" X 
which will be referred to as the polarization coordinate. We emphasize that the constructs 
V,P and T depend both on the dynamics V b and the polarization function X. They can thus 
also be viewed as images of maps V = V[Pb, X], P = P[Vb, X] and T = T[V b , X]. 



7 Function / will be referred to as middle-valued at x if f(x) = \im$_> + (f(x + S) + f(x — S))/2. Note 
that if / is continuous at x then it is also middle-valued. At points of discontinuity with finite jumps, 
middle-valued function takes a value in the middle of the jump. 
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By construction, V(X, p) carries the most detailed information about the dynamics 
Vb(qi,q 2 ) with respect to polarization properties defined by X, and will be referred to as 
the polarization representation. The other two objects are related to it via the following 
"descent" relations 

/oo pi 
d P v(x,p) = p(x) — > J dx\x\p(x) = r (15) 

i.e. P is a marginal distribution of V and T is a moment of P (and of V). In P(X) the 
information on dynamics unrelated to polarization is integrated out which corresponds to 
the X -distribution of Ref. [T] generalized to arbitrary X. X-distributions are central objects 
of interest in this work. The role of the averaged polarization characteristic < T < 1 will 
become clear in the upcoming sections. 

If X(x) G X is a one to one map of [—1, 1] onto [—1, 1], thus entailing a faithful polarization 
coordinate on the sample space, then the dynamics specified by V(X, p) is fully equivalent to 
original Vb(qi, q 2 ). The corresponding subset X-^ C X is comprised of odd, continuous, strictly 
increasing functions X(x) on [—1, 1] with unit maximal value. Moreover, utilizing such 
equivalent parametrizations is particularly efficient if they have a different iable relationship 
to base coordinates (q\, q 2 ), since then the standard analytic manipulations in the underlying 
integrals can be performed. Thus, we also consider the subset X-^ d containing those elements 
of X^ that are different iable on their domain, i.e. X^ d C C X. 

Given the utility of polarization coordinates for our purposes, it is convenient to fix 
a reference polarization representation of dynamics Vb{qi, q 2 ). Since x = X(x) G X^ d , the 
rescaled polar coordinate (|2]) is particularly suited to play this role@ Changing variables from 
(51,92) to (x,p), the associated reference polarization representation V r (x,p) is explicitly 
specified by 

/OO /"OO /»1 /»oo . 

dqi J dq 2 V b (qi,q 2 ) = J dx J dp - pV b [qi(x, p), q 2 (x, p)J (16) 



V r (x,p) 

with qi(x,p) = p cos {\{x + 1)), q 2 (x,p) = p sin (j(x + 1)). In terms of this representation, 
equations (I12H14I) simplify to 

V(X,p) = J dxV r (x,p) S(x-X{x)j (17) 
P(X) = j dxP r {x) 6(x-X{xj) (18) 



r = / dxP r (x) \ X(x) I (19) 



1 

00 



where the restricted dynamics P r (x) = / °° dpV r (x, p) focuses on polarization only, and will 
thus be referred to as the reference polarization dynamics. Note that, at the same time, P r 
can be thought of as the reference X-distribution. 



8 There are two technical reasons to choose x. First, the change of variables from base pair (91,92) to 
polarization pair (x, p) doesn't involve ^-dependence of the Jacobian. Secondly, the specification of the set 
of polarization functions in terms of x is simple. 
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To summarize, we think of all dynamics as being parametrized via distributions Vb{qi, Q2) 
in base coordinates. However, when interested in characterizing polarization properties of Vb, 
it is convenient to parametrize the same set of dynamics in the coordinate system that directly 
reflects the polarization of the sample, which is the role of V r (x,p). Integrating out the 
coordinate orthogonal to reference polarization coordinate x, we obtain a restricted dynamics 
P r (x) describing polarization only, and we thus view P r as parametrizing possibilities for 
polarization dynamics. While it is tempting to take P r (x) to be a detailed polarization 
characteristic assigned to dynamics Vb, the choice of a; as a reference polarization coordinate 
is clearly arbitrary, and we could reparametrize using any polarization function X(x) G X. 
Had we done so, V r {x,p) would have been replaced by another polarization representation 
V(X,p), P r (x) would have been replaced by X-distribution P(X), and the reference value 
of T r would have changed to T. Thus, unlike V r , P r and T r which are assigned to a given 
dynamics, we view V, P and T as "floating objects" reflecting dependence on X as well. 

A. 3 Polarization Dynamics 

As emphasized before, the objects of our interest are various X-distributions P that we can 
assign to given dynamics Vb, since they carry the information on polarization properties. 
Given that the polarization dynamics P r associated with Vb is itself an X-distribution, it is 
not surprising that, as seen from Eq. ffl8|) . the full information contained in Vb{qi, (72) is not 
needed to study its X-distributions. Rather, the polarization dynamics P r {x) represents all 
we need. We can thus write P = P[P r , X] instead of P = P[Vb, X]. 

Therefore, our analysis in what follows will focus on studying interrelations among the 
members of the triple (P r , X, P). To do that, it is desirable to be more specific about defining 
the set of possibilities for polarization dynamics. In particular, we would like to represent 
them in terms of a function set so that their description is on par with the representation 
of polarization functions via set X. In that regard, it is convenient to use the information 
stored in P r {x) in the form of the associated cumulative probability function 



which assigns probability to events [— l,x\. We will then consider polarization dynamics 
represented by S r G P where the function set P is defined as follows. 

Definition 3 (Set F ): Set ¥ is a subset of real-valued functions f\x) on [—1, 1] satisfying the 
following conditions 
(1) f(x) is non-decreasing 



(4) f(x) + f(-x) = 1 + A(x) x > -1 , A(x) = f(x) - lim^o f(x - 5) > 



Condition (4) above expresses the symmetry of the underlying dynamics in base variables. 
The function A(x) detects "atoms" in probability space associated with f(x). In particular, 
it represents the probability of a point-like event {x}, which is zero everywhere except 
possibly on a countable set of points (atoms) where it is positive. 




(20) 



(2) 
(3) 
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It should be noted that the description in terms of S r is in fact more general than the 
description in terms of P r in the sense that the density representation doesn't exist for all 
possible cumulative probability functions. Probabilistic spaces with atoms have a discrete 
part that can be described via distributions (5-functions), but even in the case when S r 
is continuous, it can be non-differentiable on a set of (Lebesgue) measure zero@ In some 
of such cases a complete description of probability space in terms of density simply doesn't 
exist. Nevertheless, we will interchangeably use both languages since we are not aware of 
a situation where this distinction would be relevant physically. More importantly, all the 
constructs considered here can be defined solely in terms of S r for all elements of P. lj Thus 
the notion of density doesn't even have to be invoked in principle. However, its convenience 
and physical utility certainly warrant its use. 

A. 3.1 Polarization Functions — Polarization Dynamics Correspondence 

For considerations that will follow, it will be useful to single out polarization dynamics 
for which S r (x) is a one to one map of [—1,1] onto [0,1]. Such subset C P consists 
simply of continuous, strictly increasing functions with unit maximal value and satisfying 
S r (—x) + S r (x) = 1. Thus, a simple transformation 

S r (x) — ► X r (x) = 2S r (x) - 1 (21) 

makes it into odd, continuous, strictly increasing function on [—1, 1] with unit maximal value, 
namely an element of X-> . The construction obviously works also in the reverse direction 
and (l2Tj) defines a natural correspondence between polarization dynamics and polarization 
functions - a bijection between P^ and 

Polarization dynamics from W do not contain atoms or events [x\,X2 > x-\\ with zero 
assigned probability. The latter can be relaxed while the above correspondence still holds 
unchanged. Formally, it is then a bijection between larger sets X c (P c ) of continuous polar- 
ization functions (dynamics). However, the extension to cases with atoms requires a minor 
modification. In particular, the problem with ( J2TT) then is that the image X r (x) is not odd 
at points of discontinuity due to property (3) of S r (x), and thus not an element of X. This 
is rectified by replacing the map (12 lj) with 

r 2S r (x) - 1 - A(x) for |ar| < 1 
S r (x) — > X r (x) = I -1 for x = -1 (22) 

I +1 for x = 1 

which only differs from ( I2T1) for S r that are not elements of P c . Map ( 1221) in fact defines a 
bijection between X and P, and we can conclude that the freedom in choice of the polarization 
functions is as large as possibilities for polarization dynamics itself. 

A. 4 Arbitrariness of Polarization Characteristics 

Viewing the constructs V, P and V as tools for characterizing the polarization properties of 
dynamics Vb with decreasing level of detail, one has to ask how meaningful these character- 
continuous functions that are not monotonic can be non-diffcrcntiable everywhere as exemplified by the 
famous Weierstrass's function. Monotonicity however greatly reduces options for such exotic behavior. 
10 This includes moments, e.g. T, as follows from properties of Sticltjcs's extension to the Riemann integral. 
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Figure 11: Possible X-distributions associated with fixed dynamics can have qualitatively 
different behavior. In this example from lattice QCD we selected two polarization functions 
from family X c (left), and two from family X^ (right). 



istics can be since they depend on the choice of the polarization function. In this subsection 
we will clarify the extent of this arbitrariness in case of X-distributions. Specifically, fixing 
the polarization dynamics via P r {x) (or S r (x)), we ask how much can X-distributions -P(X) 
(or S(X)) change as we scan the set X of polarization functions. 

To start with a relevant example, we show in Fig. [TT] several X-distributions associated 
with dynamics induced by pure-glue lattice QCD ensemble E± in two pairs of eigenmodes 
with lowest non-zero eigenvalues (see Table [1] in Sec. [3]) . The left-right decomposition of 
Dirac bispinor defines the polarization components, and the details of related calculations 
are given in Sec. [3j In the left panel we selected two polarization functions from the family 
X c of Eq. ([9]), with values of a as shown. One can see that the qualitative behavior of 
the X-distribution differs dramatically for the two choices with a = 1 case suggesting 
a double-peaked (highly polarized) behavior while a = 1/4 case exhibiting the opposite 
(highly unpolarized) tendencies. On the right panel we show a similar example for functions 
X R specified by Eq. (fTOj). 

To understand the variability of X-distributions, it is easiest to vary polarization func- 
tions within X(x) G X^* d , i.e. to consider the subset of differentiable faithful polarization 
coordinates. Then we can make the definition f|T8|) of -P(X) more explicit via analytic 
change of variable to obtain 

P(X) = — - — -p r (x-\x)) XeX fd (23) 

XMX-!(X)J V J 

where X'(x) = dX(x)/dx. The form of this result explains why we can obtain wildly varied 
X-distributions for P r (x) by changing the polarization function X(x). 

For more general analysis, it is convenient to express the underlying relations in terms of 
cumulative probability functions. In that case the defining condition for X-distribution S(X) 
associated with S r (x) and X(x) is that it is the probability of the event {x \ X(x) < X} in 
the probability space specified by S r (x). To make this more explicit, notice that if X belongs 
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to the range CR[X] C [—1, 1] of X, then {x \ X(x) < X} = [— 1, X(X)], a closed interval, while 
if X (jL Oi[X] then it is an open interval [— 1,X(X)). Here the generalized inverse X(X) of 
polarization function X(x) is defined via the relevant least upper bound, namely 

X(X) = sup{x | X(x) < X} (24) 

and it coincides with the standard inverse X _1 (X) for X G XA We then have explicitly 

{ S r (x(X)) ifXGft[X] 
S(X) = { }_ \ (25) 
I S r [X{X)-J ifXg'ftpC] 

where S r (x~) is the left limit of S r at x. The above expression assigns the X-distribution 
S(X) G P to arbitrary S r G P and X G X, and is the analog of equation f fT8]) . Taking it 
as a starting point, one can straightforwardly obtain special cases of interest. For example, 
considering faithful polarization coordinates, i.e. X G X^, we have 

S(X) = Srfx^iX)) or S r (x) = S(x(x)) XeX f (26) 

One use of the above relations is that they show, in an explicit manner, that fixing an 
invertible polarization function defines a one to one map of P onto itself, namely a bijection 
between polarization dynamics and X-distributions. Let us also restrict the polarization 
dynamics to be from P^cP (see section [A. 3. II) for which the corresponding S r (x) is a one 
to one map of [—1, 1] onto [0, 1], and thus invertible. From (|26|) it then follows that S(X) 
assigned to X G X-^ is an element of F? , and thus also invertible. We then have 

X{x) = S-^Srix)^ XeX f , S r eF f (27) 

The above equation establishes that, given a fixed polarization dynamics S r from P*, we 
can obtain an arbitrary X-distribution S(X) from F? by adjusting the polarization function 
X G X* according to (1271) . In other words, scanning through all polarization functions in X^, 
the constructed X-distributions span the whole P-' . One can continue further and make this 
result even more general, but we will not go into additional details. The conclusion is that, 
for S, G P', the arbitrariness in its X-distributions is essentially as large as the full variety 
of polarization dynamics themselves. 

The above analysis reveals the simplicity of the structure governing the interrelations in 
possible triples (S r ,X,S) when S r and X are faithful. In particular, the X-distribution S 
doesn't have to be thought of as a "dependent variable" determined by S r and X, but can be 
also viewed as an independent object from P* . In other words, we have equivalent functional 
dependencies S = S[S r ,X], S r = S r [S,X], and X = XfSV,.!?] whose explicit forms are given 
above. In what follows, we will implicitly assume the above restriction, i.e. that polarization 
dynamics are chosen from P-* , and that polarization functions are chosen from X^ so that the 
above structure readily holds. Moreover, when our discussion involves probability densities 
(P r or P), the reader can simply assume that we are working in the framework of X-^ d and 
F-f d . However, we emphasize again that everything can be straightforwardly extended to the 
P-^, X-f combination in that case. Most situations of practical interest are covered in this 
way but when the departure from P* , X-'' is necessary, we will explicitly point it out. 
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A. 5 Motivation for Relative Measures 



Considerations of the previous section suggest that, in order to meaningfully characterize 
the dynamics of polarization via X-distributions, one should proceed in one or both of the 
following ways, (i) Fix the polarization function X using well-motivated physics considera- 
tions, i.e. select polarization function that has direct physical meaning. Whether this can 
be done depends on the nature of quantity Q in question, (ii) Use only those aspects of X- 
distributions that do not depend on the choice of the polarization function. We emphasize 
that (i) and (ii) really represent conceptually different approaches to the issue. The purpose 
of this paper is to construct universal characteristics of the dynamics in the spirit of the 
latter. The following considerations turn out to be relevant for achieving that goal. 

The feature that puts all polarization functions X £ on the same footing is that 
they are equally good comparators of polarization on the sample space due to strict mono- 
tonicity. In other words, if sample (51,52) was established to be more polarized than (ti, t 2 ) 
using polarization function Xo, then this would hold also if any other X £ X^ was used to 
relate them. Thus, any characteristic that may be defined via statistics of such comparisons 
will be invariant under the choice of the polarization function. To see a simple and useful 
construct of this type, let us compare the dynamics (distribution) to dynamics vj® via 

the probability Tr that a sample drawn from is more polarized than a sample indepen- 

(2) 

dently drawn from V h . While one can approach the determination of Tr using different 



polarization functions, the answer will always be the same. Indeed, this invariance becomes 
obvious if one pictures evaluation of T R by producing simultaneous draws from and 
V\ , and recording the double-sequence of corresponding reference polarization coordinates 
{{x^\xf^),i = 1,2,.. .}. The relative frequency of encountering \xf~^\ > \xf^\ in this process 
determines Tr, and the result would not change had we tested whether |X(x^)| > |X(^- 2 ' ) )| 
instead. Obviously, the information stored in Pr l \ P^ is sufficient to carry out the above 
process and in terms of the associated X-distributions we have formally 

r R (P«,P r ( 2 ),X) = / dXdYP(X;P^,X)P(Y;P^,X) = r R (P«, Pj 2 )) (28) 

J\Y\<\X\ 

where we have made the dependence of various objects on the dynamics and the polarization 
function explicit everywhere. For obvious reasons, Tr will be referred to as the probability 
of larger polarization (PLP) in relative to . 

More detailed invariant objects can be defined using analogous reasoning. One example 
relevant for our purposes involves the probability Sr(s) that a sample drawn from is more 
polarized in the first direction than fraction 1 — s(0<s<l)of the population drawn from 
Vj® . This is invariant with respect to the choice of the polarization function because we can 
compute Sr(s) through the following process. Selecting a reference polarization coordinate 
as our measure for polarization, we again generate a double-sequence of draws from and 
Vj® . From the sequence of length N of such draws, namely {{x^i\xf^),i = 1,2,... N}, we 
can obtain an estimate Sr(s, N) of Sr(s) by first sorting {x\ } in ascending order obtaining 
a sequence {2^)} of monotonically increasing samples, and setting xq = x^Ly where i$ is 

the smallest index such that io/N > s. Note that xq approximates the separation point 

(2) 

where the population from V\ splits into fraction s of samples most polarized in the first 
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direction, and the rest. The estimate Sr(s,N) is then obtained by counting the relative 
frequency of occurrences with xf^ < Xq. It is clear that we could have used arbitrary 
X(x) G X-f to make all the above comparisons with the same result. Thus, both Sr(s,N) 
and Sr(s) = limjv-Kx) Sr(s, N) are invariant under the choice of the polarization function 
for all s. To make the contact with X-distributions, we can change the variable and transit 
from Sr(s) to Sr(X) = Sr(1/2 + X/2). This resulting object again only depends on the 
polarization dynamics involved, i.e. we have Sr(X; P^, P^). We will refer to it as relative 
X -distribution of vj, relative to 7-j . Note that, by construction, if Vf 1 = vf \ i.e. when 
comparing given dynamics to itself, we have S(s) = s and Sr(X) = 1/2 + X/2. 

A. 6 Relative Measures and their Reparametrization Invariance 

Discussion of the previous section fully defines two new objects, namely the probability 
of larger polarization Tr and the relative X-distribution Sr, both assigned to a pair of 
polarization dynamics from W. However, in order to connect them to previously defined T 
and S, we need to convert the constructive definitions into explicit ones. Let us start by a 
closer look at equation ( |28l) representing Tr which can be written as 

rVP^Pj 2 )) = f dxPW(x) ^ dyP?\y) = f dxP?\x) \ 2Sf\x) - 1 1 = 

J-l J-\x\ J-l 

= j 1 dxPj: 1 \x)\X^(x)\ = r(P«,X( 2 )) (29) 

where X^ is the polarization function assigned to dynamics (or Pr^) via natural map 
f[2T|) . Thus the PLP in P r (1) relative to P r (2) is the average polarization V of P^ in polarization 
function Xr , as defined by Eq. (fl9l) . 

To express the relations of the above type, it is convenient to use the "hat notation" for 
functional prescriptions involving domains of maps (functions) themselves. For example, the 
functional prescriptions for maps S : P x P and X : x X', introduced via 

relations fl26l) and (J2TJ) respectively, are given by 

S[S, X] = 5oI _1 and X [ S(T) , Sp) ] = S^j ° -S'(i) (30) 

where the arguments are to be thought of merely as generic elements of the function sets 
involved. It is further useful to express the natural map (I2T!) as 

5 = 5*[I]=S*oI where S*(X) = (31) 

Note that 5** G P* represents polarization dynamics associated with the uniform distribution 
P*(X) = 1/2, and that S r = S*[X r ]. According to Eq. ( 1291) we then have, in terms of 
cumulative probability distributions, that Tr = f r[ Sr , ] with 

rji[S , (i),S , (2)] = f [5(1), S^ 1 o 5(2)] (32) 
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Similarly, one can inspect that the construction of the relative X-distribution Sr corresponds 
to setting Sr = SrIS^, ] with 

Sr[S(i),S(2)] = S[S { i), S+ 1 o 5(2) ] = S {1) o o S± (33) 

Notice that since S[ S r , X r } = S[ S r , 5" 1 o S r ] = 5*, the relative construct Sr[ 5(i), 5(2) ] rep- 
resents a familiar X-distribution for dynamics 5m with respect to the polarization function 
for which the X-distribution of 5(2) is uniform. 

An important feature of the above definitions, discussed in the previous section, is that 
they are invariant under the choice of the polarization coordinate. In other words, had we 
chosen a polarization coordinate associated with polarization function X(x) instead of our 
reference polarization coordinate x, then the polarization dynamics described by Sr would 
have been described by its X-distribution S[Sr \ X] = Sr °X _1 instead. As is obvious from 
definition ( |33i) . the relative X-distribution Sr is invariant under such reparametrization. 
Consequently, being a moment of associated Pr, PLP is invariant as well. 

The basic property of PLP is that 

f r[S(i), 5(2) ] = 1 - f jR [5( 2 ), 5(i) ] (34) 

as is obvious from its statistical interpretation described in Sec. IA.51 and can be easily derived 
analytically following the definition given above. Similarly, one can immediately see from 
(1331) that the relative X-distribution satisfies 

Sr[S(2),S{i)] = S*° (Sr[S(i),S(2)]) oS 1 * (35) 

Thus, up to some necessary modifications by the 5* operation, S R [ 5(2), 5(i) ] and Sr[ 5(i), 5( 2 ) ] 
are inverses of one another, as one can easily understand from the constructive definition of 
relative X-distribution. We emphasize again that the above relations hold for Sm G W* . 



A. 6.1 Implementation 

Since the definition of absolute X-distribution, to be discussed in the next section, is based 
on the relative X-distribution, it is worth commenting on the implementation of the latter. 
The procedure is straightforward if the polarization dynamics 5(1), 5(2) (or P(\), P(p.)) are 
explicitly given. Indeed, in that case one simply implements the evaluation of 

S R {X) = 5(i)(5 ( -;(l/2 + X/2)) (36) 

or in terms of probability density 



l P i Js^ (1/2 + X/2] 

Pr(X) = f (37) 

2 P (2) (5 ( - 3 1 (l/2 + X/2; 



However, the explicit forms of relevant distributions are frequently not known in practice. 
Rather, the dynamics involved are only specified implicitly, e.g. in terms of some reduction 
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process applied to underlying more complex degrees of freedom. In such situations one 
typically has an access to the dynamics through a finite probabilistic chain of samples known 
to be generated from the associated distribution. One option to implement the relative X- 
distribution in that case is to use the constructive definition given in Sec. IA.5I An equivalent 
approach is to simply follow the original "histogram method" of Ref. [T], and apply it to the 
probabilistic chain for and the polarization function 

X(x)=2 J dyP (2) (y)-l (38) 

The specifics of the implementation then boil down to the selection of a procedure for ap- 
proximating X(x) from a finite probabilistic chain for P( 2 ). 



A. 7 Absolute X— Distribution 

The introduction of relative X-distribution doesn't change the arbitrariness of it as a tool 
characterizing polarization properties of given dynamics. Indeed, it just moves the arbitrari- 
ness in the choice of polarization function into the arbitrariness of the dynamics serving as 
a comparative standard. However, it does introduce a different angle into the problem that 
turns out to be useful. Indeed, contrary to the question "Which polarization function should 
be used?" the question "Which dynamics should we compare to?" does evoke a meaningful 
answer. 

The main idea is as follows. When seeking a characteristic describing a dynamical ten- 
dency for polarization, what we have in mind is the information on how the projection 
component behaves in relation to the complementary one, rather than what does it do "in 
isolation", or independently of it. In other words, we wish to compare to the dynamics that 
produces identical answers for any questions involving components "in isolation", but no 
effects attributable to the interplay between them. For given dynamics Vb, such dynamics 
VI is uniquely specified, and represents the behavior of statistically independent components 
with the same distributions of components themselves (marginal distributions). Specifically, 
denoting the marginal distribution associated with "P& ((71,(72) as p(q)i i- e - 

p(q) = / dq 2 V b (q,q 2 ) = / dqiV b (qi,q) (39) 
Jo Jo 

the corresponding "uncorrected distribution" V£ is given by 

K(qi,q 2 ) = p(<Zi)p(g 2 ) (4°) 

Note that the superscript "it" will be used also in other related constructs, i.e. the polarization 
dynamics associated with V^ will be denoted as P r u (or We then define the absolute 
X -distribution associated wit dynamics Vb (in its cumulative form) as 

S A = S A [V b ] = S R [S r ,S?] = SrotS^oS* (41) 

namely as X-distribution of S r relative to S 1 ". There are several remarks regarding this 
definition that might be useful to emphasize. 
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(i) The feature that makes the above construction qualitatively different from the old X— 
distribution approach is that it discards the notion of preferred fixed polarization function. 
This is obvious from definition (I4T1) corresponding to the choice of the polarization function 
X(x) = 2Sr(x) — 1, which depends on the dynamics Vb itself. More precisely, the associated 
polarization function will be the same for members of the equivalence class consisting of all 
dynamics with fixed marginal distribution in base variables. The members of the same class 
will be distinguished by their dynamical tendencies for polarization of course, namely by 
their absolute X-distributions. 

(ii) Note that the information in polarization dynamics S r is not sufficient for constructing 
the absolute X-distribution for dynamics Vb- Indeed, the marginal distribution p(q) (and 
thus S 1 ") cannot be reconstructed from S r alone. Thus, we cannot scale down from SU^ft] to 
5U[S r ]. This is related to point (i) in that for the old "fixed polarization function" approach 
this is always possible. 

(Hi) While the notion of absolute polarization function is not viable for characterizing dy- 
namical tendency for polarization, the notion of absolute X-distribution remains. Indeed, 
the construction (j4ip is invariant with respect to the choice of the polarization function to 
implement it. In other words, we could have started with other polarization representation 
leading to X-distributions S and S u for Vb, but we would nevertheless end up with the same 
absolute X-distribution Sa- At the technical level, this follows from reparametrization in- 
variance of relative X-distributions, and is made obvious by the discussion in subsection IA.5I 

A. 8 Correlation Coefficient of Polarization 

The absolute PLP is defined in the same way as the absolute X-distribution, namely 

T A = t A [Vb] = t R [S r ,S»] (42) 

and can be computed as the moment of the latter (see Eq. (Tl5l) ). Since Ta is an average 
characteristic based on comparison to statistical independence, it can be thought of as a 
correlation measure specifically designed for polarization. To be used as such, we define the 
correlation coefficient of polarization (CCP) via 

C A = 2T A -1 = S:\Ta) (43) 

Note that Ca £ [—1, 1]. The positive correlation (r^ > 1/2) signals that dynamics is acting 
to enhance polarization, while negative correlation (r^ < 1/2) indicates that dynamics works 
to suppress it. Due to its invariant nature, it quantifies an inherent property of Vb, namely 
its overall dynamical tendency for polarization. 

A natural question to ask in this regard is why the standard Pearson correlation coefficient 
is not suitable for this purpose. This is easy to see. Consider the dynamics Vb(qi, (fe) with 
support only on line q\ = q%. While the Pearson correlation will be perfect (r = 1) in this case, 
it is obvious that there is no dynamical tendency for polarization in such dynamics. In fact, 
it is just the opposite with dynamics inhibiting any possibility of polarized samples. Indeed, 
the absolute PLP is obviously zero in this situation for all dynamics with non-singular 
marginal distributions. These dynamics are thus perfectly antipolarized and characterized 
by Ca = — 1 as appropriate. 
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A. 9 Singular Dynamics 



At this point we wish to briefly comment on including in the formalism of absolute X- 
distributions the cases involving dynamics Vb for which the associated polarization dynamics 
S r is not an element of W . Let us recall that there were no restrictions on the notion of 
standard X-distribution S(X) which can be uniquely assigned to arbitrary S r G P and 
X G X. However, it was natural to focus on the case of bijective maps (elements of P^, 
X^), especially as we transitioned to the notion of relative X-distributions. This is reflected 
in our formulas for £7r[£7(i), Sq) ], such as ( 133]) and ( 1361) . by the presence of SZy In fact, 
everything discussed in the context of relative X-distributions applies to arbitrary G P 
and £7 {2) G F f . While this covers most cases of practical interest, it does not afford the 
construction of absolute X-distributions for Vb with S r £ P^ since in that case £7" is typically 
not an element of P* either. 

The main point about the extension from £7(2) G P^ to £?( 2 ) G P in the context of relative 
X-distributions is that it may not simultaneously retain the features coexisting in the former 
case. In particular, with faithful polarization comparators £7(2) G W we have the interpreta- 
tion of £7r[£7(i), £7(2) ] both in terms of standard X-distribution with particular polarization 
function, and in terms of statistical comparisons between the two dynamics, as described 
in constructive definition of Sec. IA.5I In fact, one can easily see that this dual view of the 
concept extends to the case when £7( 2 ) G P c . Indeed, the object £7~ 1 o £7( 2 ) needed in the 
definition Sr[S^, £7( 2 ) ] = £7[ £7 (1 ), £7~ 1 o £7( 2 ) ] is a valid polarization function in this case, as 
discussed in Sec. IA.3.11 At the same time, this definition is equivalent to the comparative 
one in unchanged form. 

However, the situation is different when £7(2) is discontinuous and thus represents dynam- 
ics with probabilistic atoms. The extensions in that case are non-unique in both of the two 
approaches that become non-equivalent in general. If one follows the rationale of polariza- 
tion functions, then it is natural to simply extend the S+ 1 operation using Eq. (1221) and thus 
define Sr via canonical polarization dynamics - polarization function bijection of Sec. IA.3.11 
When utilizing this definition, then absolute X-distribution assigned to dynamics that one 
generally thinks of as maximally polarized, i.e. those with P r>p (x) = \5{x + 1) + \8{x — 1), is 
P^(X) = P rp {X) and the correlation coefficient is Ca = 1- However, if one prefers to focus 
on the concept of statistical comparisons, then one is led to an alternative path which we 
now briefly describe. 

The root cause of ambiguity in the comparative definitions of Sec. IA.5I is that for distri- 
butions with atoms, comparisons using ">" and ">" can become non-equivalent, and one 
has to fix the treatment of the "=" case. To see the associated issues, let us again consider 
the polarization dynamics P r>p , but assume in addition that the underlying dynamics Vb 
has a marginal distribution p(q) = %5(q) + \5{q — 1). This implies that its uncorrelated 
polarization dynamics is P™ p (x) = \5{x + 1) + \8{x) + \8{x — 1). Adopting the compara- 
tive scheme with sharp inequalities (as we did in Sec. IA.5j) . one immediately obtains that 
^r[Pt,pi Pr, P ] = 1/2 and Tr[P^ p , P np ] = 0. On the other hand, admitting the equality always 
as a positive outcome, we obtain f R[P r , p , = 1 and Tr[P^ p , P T)P \ = 1/2 instead. Note that 
in both cases the equation reflecting conservation of probability is not satisfied. In the 
former case, some probability is "lost" by ignoring the equal outcomes, while in the latter 
case it is "created" via double counting. The natural approach is to insist that Eq. (IM)) 
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Figure 12: Construction of absolute X-distribution for singular dynamics characterized by 
Pr,p(%) ( see text) in the context of statistical comparisons. The intermediate step, Sr(s), is 
shown on the left while Pa(X) is on the right. 

be satisfied and that the two dynamics being compared are treated on equal footing. This 
forces us to assign half of the equal outcomes in favor of first dynamics and the other half 
to the second dynamics. That way we obtain T R[P r , p , P" p ] = 3/4 and Tr[P^ , P r>p ] = 1/4, 
implying that Ca = 1/2 in this case. We emphasize that while we are focusing on a specific 
example, the above rule uniquely fixes the assignment of PLP and the correlation coefficient 
for all possible dynamics. 

The above approach can be consistently extended to the level of X-distributions as 
well. As one can easily check, the construction of Sec. IA.5l for relative X-distributions only 
becomes ambiguous when both dynamics have an atom at the same point Xq. The proper 
prescription is that the segment of cumulative probability function associated with such 
points be linear. Illustrating this again on the case of P r<p and P" we show in Fig. [12] the 
construction of Sr for P rjP relative to P" and the associated absolute X-distribution P4. 
Note that J dX Pa(X) |X| = 3/4 = Ta consistently with the above prescription for PLP. 

A. 10 More on the Invariance of Absolute Measures 

Apart from reparametrization invariance, absolute measures possess additional symmetry 
that we would like to mention. To do that, it is instructive to think of reparametrization 
invariance in the following way. Imagine that the generator based on the dynamics Vb{qi, Q2) 
produces the sequence of samples {(<7i, ^2)1}, distributed accordingly. The operator of the 
device broadcasts the sequence to various distant receivers, each of which has a polarization 
meter at his disposal. When fed a sample (qx, q 2 )i, the polarization meter outputs its degree 
of polarization Xj, thus allowing each receiver to analyze the polarization properties of P&. 
While the operation of each existing polarization meter is based on a valid polarization 
function 1 6 X^, there is no standard choice adopted by the manufacturers. Thus, each 
receiver obtains a different sequence {Xj} to work on in general. This would seem to preclude 
the possibility of producing equivalent polarization characteristics across the spectrum of 
receivers. 
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However, every receiver is required to follow special instructions, identical for all, on how 
to use his polarization meter to determine the dynamical polarization properties of the trans- 
mitted dynamics. According to this recipe, in addition to running {(qi, 52)4} through the po- 
larization meter and producing {Xj}, he is also required to form an independent sequence lj 
{(?i> ^2)^}) i n which the received components are assigned to one another randomly rather 
than as dictated by the dynamics. The receiver is then to run this randomized sequence 
through his polarization meter to obtain the associated {X^}. Finally, he is instructed to 
perform the statistics of comparisons between and Xf to yield the desired dynamical 
measure of polarization (absolute X-distribution or correlation coefficient). Reparametriza- 
tion invariance expresses the fact that, assuming they follow the above recipe, all receivers 
will end up with exactly the same dynamical polarization characteristic. 

Thus, by virtue of providing the above instructions, the sender doesn't have to worry 
about the fact that receivers use various kinds of polarization meters. The additional sym- 
metry of absolute measures pertains to the fact that, to some extent, the same is also true 
in reverse. In particular, by virtue of the same prescription, the receivers are protected from 
the fact that the sender could have used various kinds of generators. Indeed, recall that the 
purpose of the formalism that we developed was to extract from the full dynamics of quan- 
tity Q the information on polarization. In doing that, we postulated that the information 
on dynamics of magnitudes qi, qi of the projections, i.e. Vb(qi, 92), is sufficient for this since 
the relevant measures are based on the ratios of these magnitudes. However, we could have 
used the distribution of (qf , q%) (a > 0) associated with the same dynamics of Q for this 
purpose. Thus, different senders communicating information on the same dynamics could 
have used different kinds of generators and sending the sequence of samples {(qf, q%)i}, in- 
stead of {(q±, 92)1}- As one can easily see though, the above recipe for constructing absolute 
measures guarantees that the result is insensitive to this freedom because the statistics of 
required comparisons will not change. Formally, these measures are invariant under 



It is due to simultaneous validity of both types of above invariances that we view the con- 
structs discussed here as universal polarization characteristics truly characterizing the dy- 
namics itself, rather than representing our choices for description of polarization. 

B Elements of Implementation 

In this Appendix we will discuss certain numerical aspects of the results presented in Sec. EJ 
While the implementation described below was applied specifically to the case of overlap 
Dirac eigenmodes, one can of course use these strategies in more general settings. We will 
pay most attention to the computation of absolute X-distributions, but will also discuss the 
calculation of correlation coefficient Ca and the procedure for the determination of chiral 
polarization scale At- 

Before we start, it is useful to recall that the input for our calculations is certain dynamics 
^6(91,92)) but we only have a restricted statistical knowledge of it from relevant numerical 

11 What we mean by "independent" in this case is independent with respect to the ordering of the elements. 
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simulation. In other words, what we have is a finite population of pairs { (q^ , q\ ) | % = 
l,...,N}, whose statistics is governed by Vb(qi, q 2 )- By virtue of this finite population, the 
latter can be thought of as being "approximated" by a discrete distribution 

1 N 

V b (qi,q 2 ) — ► V b (q 1 ,q 2 ) = —^26(q 1 -q[ l) )5(q 2 -q2 ) ), 

i=i 

N N 

V^q 2 ) — ► ^(ft.^^-^^-gW)-^^-^)) (45) 

4=1 i=l 

where the second equation relates to the "uncorrected" counterpart V^iqi, q 2 ) of Vb(qi, q 2 )- 
In case of Dirac eigenmodes the population comes from local values in the selected group 
of modes, such as those at scale A in given gauge ensemble, with pairs representing the 
magnitudes of right and left chiral components. 

Our implementation followed the primer calculation described in Sec. [2J Technical details 
not specified there have to do with the finiteness of the population representing Vb(qi,q 2 ) 
and the fact that, while one wishes to have its size N as large as possible in order to represent 
the dynamics adequately, the computational complexity involved in handling the complete 
information as expressed in Eq. ( I4"5j) grows as N 2 in case of uncorrelated distribution. Since 
N can be quite large even with relatively modest size of gauge ensembles, i.e. we have 
iV w 6 x 10 7 for ensemble E± with only two eigenmodes per configuration contributing, it is 
computationally advantageous to do some coarse-graining as explained below. 

Absolute X— distribution. Given our statistics, we have decided to evaluate Pa(X) 
at Nfy = 50 equidistant values based on the information stored in the associated AT& = 50 
bins of probability. This allows us to plot the distribution with enough resolution to capture 
relevant details while it also allows for enough statistics in each bin so that the statistical 
errors are reasonably small. In order to reduce the computational effort, we organized the 
calculation as follows: 



1. Compute the reference polarization coordinates xi, x^-i representing the bound- 
aries of each bin. To avoid the numerical evaluation of inverse tangent and using the 
fact that our construction doesn't depend on the parametrization of the polar angle, 
we actually work directly with slopes t±, t^-i (t = q 2 /qi)- These slopes separate 
the sample space quadrant into segments containing equal number of samples drawn 
from the uncorrelated distribution V^. In Fig. [1] we indicated with solid gray lines 
these boundaries when the quadrant is partitioned into 10 such bins. 

2. Count how many samples (q[ , q^ ) drawn from the correlated distribution Vb fall into 
each bin by comparing their slope with the boundaries determined in the previous step. 

Since we only need to inspect each sample once to determine which bin it belongs to, 
the second step has a complexity O(N), which is negligible compared to the first step which 
is 0(N 2 ) in the naive implementation. One possibility to deal with the first step is to use 
a randomly drawn sub-population from the corresponding N 2 samples to approximate the 
uncorrelated distribution. The problem with this approach is that it interferes with the 
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error analysis. However, one can organize the calculation carefully and compute the first 
step using a method that has a complexity only O(NlogN). 

The key to this is the ability to compute S^(t) efficiently Indeed, the problem of finding 
the boundary is equivalent to solving the equation S^(ti) = i/Nb, and one can solve this 
equation using a bisection method which can determine U to one part in a million with 
only 20 evaluations of S^(t). To compute S^(t) we sort {q^ } in ascending order, a task 
of complexity O(NlogN). For every we then compute the fraction f\ of elements from 
{ } smaller than tq[ l \ Since {q^ } is sorted, one can use a binary search which has a 
O(logX) complexity. Computing the fraction /j for all X values of to get the answer 
S™(x) = fii the total complexity for this task is 0(N log N) as advertised. 

While the above strategy makes the calculation of X-distributions manageable, it can still 
take a considerable amount of time. To speed it up, we sort the values q^ and separate them 
into bins with equal number of samples. Each bin is then replaced by a single representative 
value equal to the average of the samples it contains. For the smallest ensemble E 1 where one 
can comfortably carry out the calculation without binning, we compared the results of the 
binned version using 1000 elements per bin with the exact calculation. We found that the 
differences are much smaller than the error bars. For larger ensembles, the approximation 
is expected to be even better since the number of samples is increased. The binned version 
was thus used in our evaluations. 

The result of the above calculation determines the fraction of samples drawn from Vb 
belonging to each of the Nf, bins. This information allows us to plot the distribution (his- 
togram) of dS r /dSr as a function of S*", which can then easily be scaled into the absolute 
X-distribution Pa(X). To estimate the errors for each value of argument X (or bin), we 
used single elimination jackknife method with new samples generated by removing one con- 
figuration from the gauge ensemble. 

The Correlation Coefficient. The evaluation of correlation coefficient Ca is equivalent 
to the calculation of Fa, namely the probability of larger polarization. While r a can be 
evaluated independently of the absolute X-distribution, we follow its definition as a moment 
of Pa{X) since this streamlines the error analysis. Every jackknife sample (PaY in the 
process of calculating the absolute X-distribution produces the jackknife sample P^. In 
particular, if (Pa)1 is the value of the estimate for bin i, then we have explicitly 



where Xj are the boundaries of the bins that divide the domain [—1,1] into Nb equal size 
bins, and thus X = — 1 and X Nb = 1. The jackknife analysis of the above estimates then 
determines the mean value of the correlation coefficient and its error. 

Chiral Polarization Scale. One of the main results in this work is that the correlation 
coefficient Ca changes sign at a specific scale A T of the eigenmodes. If we denote by Ca(A), 
the correlation coefficient evaluated using a pair of eigenmodes Ai ; 2 "bracketing" the scale 
A, i.e. Ai < A < A2, then the chiral polarization transition point is defined via Ca(^-t) — 0. 
To determine this scale, we first crudely map out the behavior of CU(A) by evaluating the 
correlation coefficient on a grid with steps of about 200 MeV. Note that, using the spectral 




(46) 



i=l 
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information from Table [I] we choose the overall position of these grid points such that all 
gauge configurations in the ensemble contribute modes for the calculation. Determining 
the transition region this way, we then evaluate Ca(A) at few more points in the estimated 
vicinity of the transition point. Using this information, several values Aj closest to the 
transition are finally selected. The correlation coefficient is then evaluated for each of the 
jackknife samples and the estimate of the correlation matrix for CU(Aj), using a jackknife 
method, is formed. The final result for A^ and its error are evaluated using a linear regression 
based on this correlation matrix. 

Calculation of Eigenmodes. We determine the eigenvalues of the zero-mass overlap 
operator using our own implementation of the implicitly restarted Arnoldi algorithm with 
deflation [12]. To speedup the calculation and to reduce the amount of memory required, 
we compute the eigenvalues of D^D = + D which is Hermitian and commutes with 75. 
This allows us to project on a chiral subspace and compute only one member of the chiral 
pair. The eigenvectors and eigenvalues of the overlap operator are then straightforwardly 
constructed using the basic properties of the overlap Dirac matrix. 
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